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lHj' ABSTRACT 

We use WIRC, IR images of the Antennae (NGC 4038/4039 ) together with the exten- 
sive catalogue of 120 X-ray point sources ()Zezas et al.ll2006l) to search for counterpart 
candidates. Using our proven frame-tie technique, we find 38 X-ray sources with IR 

■ counterparts , almost doubling the number of IR counterparts to X-ray sources first 

identified in IClark et all (|2007f ). In our photometric analysis, we consider the 35 IR 
counterparts that are confirmed star clusters. We show that the clusters with X-ray 
sources tend to be brighter, K s k, 16 mag, with (J — K s ) = 1.1 mag. 

We then use archival HST images of the Antennae to search for optical counter- 
parts to the X-ray point sources. We employ our previous IR-to- X-ray frame-tie as an 
intermediary to establish a precise optical-to- X-ray frame-tie with < 0.6 arcsec rms 
positional uncertainty. Due to the high optical source density near the X-ray sources, 
we determine that we cannot reliably identify counterparts. Comparing the HST po- 
sitions to the 35 identified IR star cluster counterparts, we find optical matches for 
rS | 27 of these sources. Using Bruzual-Charlot spectral evolutionary models, we find that 

■ most clusters associated with an X-ray source are massive, ~ 10 6 Mq, young, ^10 6 

yr, with moderate metallicities, Z — 0.05. 
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1 INTRODUCTION 

The numerous X-ray point sources and young, massive star 
clusters in the Antennae make this galaxy pair an ideal 
laboratory for studying the environments of X-ray binaries 
(XRBs) . Previously, Chandra observations revealed 49 X-ray 
poin t sources ranging in luminosity from 10 38 — 10 40 ergs 
s" 1 (|Zezas et all 120021 ) . A 411 ks total exposure, consisting 
of six additional pointings spre ad over two years , revealed 
an additional 71 X-ray sources IjZezas et al. 2006) down to 
a luminosity of 2 x 10 37 ergs s _1 . While most are XRBs 
with either a black hole or neutron star, those sources with 
Lx > 10 39 ergs s _1 are more unusual objects classified as 
ultraluminous X-ray sources (ULXs). Some theories suggest 
these ULXs are intermediate mass black holes with masses 



from 100 - 10,000 M n (e.g. iFabbianol Il989l; IZezas et al.l 
1 19991 : iRoberts fc WarwicM2000l :l Makishima et al .1120001 ) but 
there remains a consid erable amount of controversy con cern- 
ing their origins (e.g., iKing et al]|200ll ; lRobertsll2007l . and 
references therein). 

For this study, we use a distance to the Antennae of 
19.3 Mpc (for H = 75 km s" 1 Mpc" 1 ), as determined 
by the spectroscopic redshift (|Zezas et al.| [2006). However, 
there exists some ambiguity concerning the distance to this 
galaxy pair. While a distance derived from the Type la 
supernova (SN) 2007s r suggests a distance of 22±3 Mpc 
l|Schweizer et al,ll2008l ). analysis of the RGB colors indicate 
that t he distance can ran ge from 22 Mpc to as low as 13.8 
Mpc l|Saviane et al.l [20041 ) . While this lower distance would 
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halve our luminosity and mass estimates, this factor of a 
few will not affect any statistical study of their population. 
Hence, we do not anticipate these discrepences in the dis- 
tance to NGC 4038/4039 as being problematic. 

Compact objects tend to be associated with massive 
star formation, which som e theories suggest is predominant 
in young stellar clusters (|Lada fc Ladal I2003T ). In previous 
work on the Antennae we find a close association between 
compact objects and clusters, identifying 15 p ossible IR 
counterparts to X-ray sources |Clark et alJl200TT ). Many of 
these counterparts reside in the spiral arms and 'bridge' re- 
gion of the Antennae - locations predominant in massive star 
formation. The spiral arms wrap around the northern nu- 
cleus and extend to the south of the southern nucleus, while 
the bridge region is the dense area connecting the two galax- 
ies. Those X-ray sources without counterparts could be com- 
pact objects that escaped their parent cluster or remained 
behind after their cluster dissolved. In IClark et al l (2007) 
we suggest a third possibility, that some X-ray sources do 
have counterparts, but these are too faint to see in the IR 
images. 

In this paper, we expand on our initial study by search- 
in g for IR counterpa rts to the 120 X-ray sources identified 
in lZezas et al.l (2006). We then extend this work to optical 
wavelengths using HST images of the Antennae. The higher 
sensitivity of HST in non-dust obscured regions allows us 
to search for additional counterparts to X-ray sources. In 
addition, combining multi-band photometry in the optical 
and IR, we can use spectral evolutionary models to measure 
cluster properties. In §2 we discuss our IR observations and 
analysis of optical HST archival images. We describe our 
photometric analysis of counterpart cluster properties in §3, 
and summarize our results in 84. 



2 OBSERVATIONS AND DATA ANALYSIS 

2.1 Infrared and Optical Imaging 

We base this study on near-infrared (IR) and optical HST 
images of the Antennae. We r eport the red u ction and anal- 
ysis of the near-IR images in iBrandl et alj (|2005l ). In sum- 
mary, we acquired 20-minute total exposures in the J (1.25 
(jm) and K s (2.15 fim) filters using the Wide-field InfraRed 
Camera (WIRC) on the Palomar 5-m telescope during the 
night of 2002 March 22. We obtained the optic al HST im- 
ages from the literature (jWhitmore et al]|l999l ). This data 
set was acquired using the WFPC2 camera and consists 
of images in the following filters: F336W((7), F439W(S), 
F555W(U), and F814W(/). 

In our efforts to understand the environments of the 
Antennae X-ray sources, we made frame-ties across three 
wavelength regimes: Chandra X-ray, WIRC IR and HST op- 
tical images. Using the IR as an intermediary between the 
optical and X-ray frame-tie has many advantages over direct 
optical/X-ray matches. First, the IR penetrates dust in sim- 
ilar ways to X-rays, facilitating the identification of counter- 
parts to X-ray sources. Furthermore, while bright IR sources 
may be obscured in the optical, the converse is generally un- 
true - any bright HST source shows up prominently in the 
IR images, enabling an excellent optical/IR frame-tie, and 
thus closing the astrometric loop at all wavelengths. Previ- 
ous attempts to match X-ray and optical positions produced 



many possible counterparts, but with poor reliability — as 
many as 75 percen t of the matches are chance coincidences 
(|Zezas et al.ll2002l ). 

In the new catalogue of X-ray sources, IZezas et alj 

(2006) revised the X-ray source positions and the new posi- 
tions differ by 0.2 ±0.4 arcsec in R.A. and 0.6 ± 1.8 arcsec in 
Dec a s compared to those coordinates listed in lZezas et alj 
(2002). Considering we are searching for IR counterparts 
within 1 arcsec of an X-ray s ource position (se e below), we 
adjusted our initial frame-tie (|Clark et al.ll2007T ) between the 
WIRC, IR j ^ 3 -ban d image and the X-ray sources listed in 
IZezas et al.1 (|2002|). to account f or the new, X-ray coordi- 
nates listed in lZezas et alj |2006l ). This modification to the 
initial frame-tie yielded the same rms positional uncertainty, 
~ 0.5 arcsec. 

Once our frame-tie between the IR and X-ray was in 
place, we extrapolated X-ray positions to optical HST I- 
band positions. Tying Chandra X-ray coordinates directly to 
HST positions is nontrivial due to field crowding in the HST 
images. Instead, we used our excellent frame-tie between the 
IR and X-ray as an intermediary. Selecting only bright, iso- 
lated optical sources, we made our frame-tie using 12 targets 
corresponding to sources seen in the IR. Using the mapping 
method describe above, we matched IR pixel positions to 
HST right ascension and declination. This frame-tie yielded 
an rms positional uncertainty of less than 0.6 arcsec. Apply- 
ing our IR-to-optical astrometric fits to previo usly derived 
ifs-b and, x,y pixel positions of X-ray sources |Clark et al.1 
l2007h . we found the HST positions for all Chandra X-ray 
sources. 



2.2 Identification of Infrared Counterparts to 
Chandra Sources 

Once our astrometric frame-tie was in place, we found a to- 
tal of 28 likely and 10 possible IR counterparts to Chandra 
X-ray sources in the Antennae, where likely counterparts 
are defined to be within a 1.0 arcsec radius (equivalent to 
twice the positional uncertainty for our frame-tie) of an X- 
ray source and possible counterparts are between 1.0 and 1.5 
arcsec (four-to-six times the positional uncertainty) from an 
X-ray source. The 35 sources that are confirmed star clusters 
are listed in Table 1 and shown in Fig. 1, with subimages 
of each source displayed in Fig. 2. As can be seen in the 
figure, the counterparts are almost exclusively in the spiral 
arms and bridge region. Of the 38 X-ray sources with coun- 
terp arts, two are the nu clei (X-50 and X-61 as listed in table 
5 of IZezas et al.ll2006l ). one is a background quasar (X-90; 
IClark et al. I l2005l ). and two lack measured X-ray luminosi- 
ties (X-24 and X-52). The source X-107 is on the edge of the 
frame, making sky-subtraction difficult, so we excluded it. 
Therefore, in our analysis of cluster properties, we only con- 
sidered the 32 IR counterparts that are confirmed star clus- 
ters and that have X-ray sources with measured luminosi- 
ties. Furthermore, of the 15 IR cluster counterparts found in 
IClark et al.l |2007l ). four drop out as counterparts in this new 
analysis — the X-ray sources X-10, X- 26, X-15, and X-22 , 
following the numeration in table 1 of IZezas et al.1 (|2002l ). 

are no longer associated with a star clu ster. 

Using the technique developed in IClark et al.1 l|2007h . 
which is based on the IR source density in the field around 
each X-ray source position, we attempted to estimate the 
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Table 1. IR Cluster Counterparts to Chandra X-Ray Sources. Act and A5 are the positional 
offsets in units of arcsec from the Chandra coordinates to the WIRC coordi nates of the p r opose d 
counterpart. Chandra Src ID numbers follow the naming convention of IZezas et al.l ||2 006), 
The listed RA and Dec Chandra coordinates have an uncertainty of 0.5 arcsec llZezas et al,l 
2006). The near-IR J and K s magnitudes were measured using PSF photometry. The values in 
parenthesis for each magnitude are uncertainties in the final listed digit. 



Chandra Src ID 


RA 


Dec 


Act (arcsec) 


A5 (arcsec) 


J 


K s 


Likely 














Counterparts 














7 


1 9-01 -4Q 6zL 


1 8-^9-07 1 n 


0.45 


0.66 


1 Q Q^ff) 91 1 


1 £ XS/Tl 1 R\ 

LO.OO^U. -LO j 


\\ 


1 9 01 -50 3^ 

. u j_ . 'jyj . o o 


-1 8-52-1 ^ 50 

lu. Ji. i-'J . ou 


0.30 


0.76 


1 fi 47(0 02"! 


1 ^ XQ(T) 09^ 


16 


1 2 01 -50 51 

Li.UL.ij'U.OL 


-1 8-52-04 30 

LO . t^ii . Urt . iJU 


0.15 


0.08 


i fi 21 mi 

LU.iLl U.UL J 


1 fififO 021 

LiJ.UUI U.Ui j 


22 


1 201 -51 01 


-1 8-52-33 00 


0.45 


0.55 


1 7 ^fifO 02*1 

l i . y u. L/ij j 


1 6 59C0 01 "1 

LU.UC/l U.UL ) 


24 


L^ . U L . O L . Trt 


- LO . 0.£ . OL . 1 U 


0.15 


0.18 


1 7 ^fifo r\9\ 

L I .0\J\\J.U^ J 


1 fi ^Q/T) 01 ~\ 
LU . ijy^u. u L y 


27 


1 901 -51 32 


lOi Oi . il U . TtU 


0.45 


0.42 


1 X 97f'0 01 ^ 

lO.i 1 1 U.UL I 


1 7 ^7(0 08 s ! 

LI . O I v U . UO J 


38 


1 9 01 -52 02 

_1_ i . U L . < J i . U i 


-1 8-52-27 50 


0.60 


0.30 


1 fi 82(0 01 "l 


1 6 1 0f0 01 "1 

LU.1UI U.U1 j 


47 


1 9-01 -^9 7^ 

L^ . \J L . Oi . I O 


-_LO.tJi.OU. LU 


0.15 


0.09 


is 4&(o r\d\ 


1 7 78(0 09~\ 

L 1 . ( O^U.U^ ^ 


51 


1 2 01 -53 00 


-1 S- r >9-0Q 1 n 


0.75 


0.04 


1 5 77(0 01 ~t 

LU. 1 I 1 U. Ul 1 


1 ^ 1 SCO Ofi 1 ) 

LiJ. 1 Jl U.UU 1 


52 


1 9-01 OS 


1 S- c i9-9'^ £0 


0.60 


0.52 


17 1 f\(C\ 09^ 


1 fi 91 ('n oi ^ 


58 


12:01:53.41 


-18:53:07.00 


0.45 


0.20 


15.79(0.18) 


15.51(0.10) 


60 


12:01:53.42 


-18:53:51.30 


0.45 


0.62 


21.95(1.68) 


17.97(0.09) 


83 


12:01:54.50 


-18:53:06.70 


0.00 


0.57 


16.71(0.08) 


16.45(0.07) 


84 


12:01:54.56 


-18:53:03.80 


0.15 


0.00 


15.05(0.03) 


14.27(0.02) 


85 


12:01:54.62 


-18:52:09.40 


0.45 


0.05 


18.52(0.20) 


16.63(0.01) 


86 


12:01:54.78 


-18:52:52.00 


0.15 


0.11 


16.76(0.02) 


14.88(0.04) 


87 


12:01:54.84 


-18:52:14.80 


0.15 


0.25 


16.60(0.03) 


15.92(0.01) 


88 


12:01:54.98 


-18:53:06.60 


0.30 


0.03 


16.52(0.01) 


14.66(0.01) 


91 


12:01:55.05 


-18:52:41.00 


0.90 


0.50 


17.81(0.01) 


16.89(0.02) 


94 


12:01:55.17 


-18:52:47.90 


0.45 


0.09 


17.10(0.07) 


15.71(0.02) 


98 


12:01:55.54 


-18:52:23.80 


0.30 


0.03 


15.86(0.07) 


15.21(0.04) 


99 


12:01:55.65 


-18:52:15.20 


0.60 


0.13 


16.30(0.04) 


15.40(0.01) 


101 


12:01:55.74 


-18:52:42.70 


0.15 


0.57 


16.70(0.02) 


16.10(0.02) 


102 


12:01:55.74 


-18:52:06.10 


0.45 


0.01 


17.81(0.02) 


17.48(0.02) 


107 


12:01:56.47 


-18:54:42.00 


0.15 


0.06 






117 


12:01:58.42 


-18:52:49.60 


0.30 


0.80 


17.64(0.02) 


16.90(0.04) 


Possible 














Counterparts 














14 


12:01:50.47 


-18:52:21.80 


0.75 


1.01 


15.86(0.03) 


15.12(0.02) 


15 


12:01:50.47 


-18:52:12.70 


1.50 


0.29 


14.98(0.02) 


14.06(0.01) 


36 


12:01:51.85 


-18:52:27.80 


0.75 


1.25 


16.77(0.04) 


15.94(0.01) 


41 


12:01:52.19 


-18:52:20.60 


1.35 


0.51 


17.06(0.04) 


16.18(0.02) 


81 


12:01:54.35 


-18:52:10.30 


1.20 


0.64 


20.30(0.45) 


16.84(0.03) 


89 


12:01:54.96 


-18:52:32.50 


0.15 


1.19 


18.11(0.02) 


17.47(0.04) 


92 


12:01:55.06 


-18:52:16.50 


0.90 


1.15 


17.25(0.03) 


15.77(0.01) 


95 


12:01:55.37 


-18:52:48.90 


0.15 


1.08 


16.21(0.02) 


15.27(0.04) 


119 


12:01:59.20 


-18:51:37.30 


0.60 


1.01 


17.78(0.01) 


15.97(0.01) 



level of 'contamination' of these counterpart samples due 
to chance superpositions of unrelated X-ray sources and IR 
clusters. We expect 8 sources with a la uncertainty of +0.1/- 
O.lQ of the 28 likely counterparts to be due to chance super- 
positions of unrelated objects, and 10 with a la uncertainty 
of +0.4/-0.3 1 for the 10 possible counterparts to be chance 
superpositions. 

In Fig. 3 we plot the X-ray luminosity versus the sep- 
aration in arcseconds to the associated IR cluster counter- 
part. We also divided the X-ray sources into three luminos- 
ity bins: Low Luminosity X-ray sources (LLX's) had Lx < 



1 Found using confi dence levels fo r small number statistics listed 
in tables 1 and 2 of lGehrelsl Jl986f) 



3 x 10 38 ergs s -1 , High Luminosity X-ray sources (HLX's) 
were between Lx of 3 x 10 38 ergs s" 1 and 1 x 10 39 ergs s -1 , 
while Lx > 1 x 10 39 ergs s _1 were Ultra-Luminous X-ray 
Sources (ULX's). A ll luminosities were taken from table 5 of 
IZezas et al.l (|2006h . Notice that there is no trend in separa- 
tion from the IR cluster counterpart and X-ray luminosity. 
This seems to indicate no preference as to where XRBs of 
different luminosity classes form in star clusters. 

2.3 Identification of Optical Counterparts to 
Chandra Sources 

The complex field seen in the HST images of the Antennae 
makes finding counterparts to X-ray sources difficult and ne- 
cessitated a different method for defining counterparts than 
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Figure 1. X-ray positions from Zcz as et al. (2006) overlaid on the WIRC, K B band image of the Antennae. The right panel is a 
1.6 X 2.0 arcmin blow-up of the central regions of the galaxy pair showing the X-ray source positions in more detail. Red circles are 
those X-ray sources with IR counterparts, while black circles are X-ray sources lacking counterparts. 



that used for the IR (see above and IClark et al.ll2007h . Us- 
ing our precise frame-tie we defined areas of positional un- 
certainty around each X-ray source. Specifically, an inner 
aperture with a radius of 1 arcsec and an annular region 
from 2.0 - 3.0 arcsec. We then defined possible matches as 
those X-ray sources with two optical sources in the circular 
aperture and less than five optical sources in the annulus, 
and likely matches as those X-ray sources with one opti- 
cal source in the circular aperture and less than five optical 
sources in the annulus. If more than five sources lay in the 
annular region, we considered the region to be too compli- 
cated for a positive counterpart identification, regardless of 
how many sources lay in the inner aperture. Using these cri- 
teria, we identified seven X-ray sources with likely matches 
to a single /-band source and one X-ray source, with possi- 
ble matches to two /-band sources. We used the HST /-band 
image to search for counterparts because this filter covers 
the longest wavelength of the available HST bands and so is 
least affected by extinction. 

Repeating the procedure discussed in IClark et al.l 

(|2007l ). we estimated the level of source contamination as- 
sociated with our identified optical counterparts to X-ray 
sources. We expect five with a la uncertainty of +0.5/-0.3 1 
of the seven likely counterparts to be due to chance superpo- 
sitions of unrelated objects, and seven with a la uncertainty 
of +3.8/-5.4 1 for the one possible counterpart to be chance 
superpositions. Clearly these statistics indicate the major- 
ity of our optical counterparts are chance superpositions and 
further demonstrate the difficulty of making such matches 
in the complex structure of the HST images of the Anten- 
nae. Therefore, we did not perform a photometric analysis 



on these source as we could not reliably identify the counter- 
parts and so could not provide any statistically meaningful 
information on the X-ray source environments. 

Instead, we considered the optical equivalents to the 32 
IR cluster counterparts identified in this work. These optical 
matches were identified using the IR-to-optical astrometric 
frame-tie to match IR counterpart positions to HST posi- 
tions. As we discuss below, in many cases a single IR coun- 
terpart split into multiple optical counterparts, and we la- 
belled these conglomerations as a positive match. We found 
optical counterparts to 27 IR cluster counterparts and Fig. 
2 displays subimages of those counterparts to X-ray sources 
seen across all six IR and optical bands. Fig. 4 is an /-band 
image of the Antennae showing the positions of all counter- 
part candidates to X-ray sources. 



3 RESULTS AND DISCUSSION 

3.1 Photometric Properties of the IR 
Counterparts 

3.1.1 Colour Magnitude Diagrams 

We made (J — K B ) versus K s colour magnitude diagrams 
(CMDs) using the 226 star clusters with J and K 3 mag- 
nitudes. T he cluster m a gnitu des for 219 sources were first 
reported in lClark et alj l|2007h and this study includes pho- 
tometry for an addition al seven clus t ers. F ig. 5 expands on 
the CMD presented in IClark et al.l l|2007h . by including a 
larger sample of cl usters associated with X-ray sources. As 
we did previously (|Clark et alj 120071 1 . we divided the X-ray 
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Figure 2. (left panel) Subimages from the K s band image of the Antennae for X-ray sources with cluster counterparts. Each subimage 
is 10 X 10 arcsec, centered on the cluster counterpart and the circles are 1 arcsec in radius, centered on the X-ray source position. X-ray 
sources labelled in red show possible matches to cluster counterparts, (right panel) Subimages highlighting optical counterparts to 
X-ray sources with IR cluster counterparts. The small circles are positional error circles with 1.0 arcsec radii centered on X-ray source 
positions. The larger circles are centered on cluster counterparts and have 1.3 arcsec radii equivalent to the i^s-band photometric 
aperture. We label those X-ray sources with possible counterparts in red. Each image is 6.0 X 6.0 arcsec. 



6 D. M. Clark, et al 



0.5 



0.0 



— 1 — 1 — 1 — 1 — 1 — 1 — 
A 


i | i 
A A 


11,111 

A 

A 


i i i i | i i i i 


A 




A 




A 

A A 
A 


A 

A A 


A 1 
A 

A 

: A 
A A i 


A A 


A 


A 


A 




A 

i i i i 1 i 


, , 1 . 


AA A A 
, . . : l , . . , 





38.0 38.5 39.0 

Log(L) 



Figure 3. Plot of X-ray luminosity versus separation from the IR 
counterpart. The dashed lines divide the plot into three separate 
regions: LLX: L x < 3 X 10 38 ergs s _1 , HLX: 3 X 10 38 ergs s _1 
< L x < 1 X 10 39 ergs s" 1 , and ULX: L x > 1 X 10 39 ergs s _1 . 
Notice that there is no trend in X-ray luminosity, nor luminosity 
class, in terms of separation from an IR counterpart. 



Table 2. Summary of Potential IR Counterpart Prop- 
are uncertainties in each value. 
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sources into the three luminosity classes, LLX, HLX and 
ULX (see §2.2 above). As can be seen in the plot, the ma- 
jority of the X-ray sources with IR counterparts are LLX's, 
26, while there are five HLX's and two ULX's. Intere stingly, 
in the new list of X-ray sources (Zez as et al.l |2006|'). three 
that a re associated with IR star clusters in both Clark et al.l 



and here ch ange in X-ray lumi nosity class from those 
published values in lZezas et al. [ ||2002D ; X-25 goes from HLX 
to LLX, X-32 from ULX to HLX, and X-40 from LLX to 
HLX. Comparing the old and new sample of IR counter- 
parts, the new sample is about a half a magnitude fainter 
(see Table 2), but shows no difference in colour compared 
with our previous study (|Clark et al.ll2007l table 4). 



3.1.2 Cluster Mass Estimates and r\ 

Following the procedure outlined in IClark et al.l (|2008l ). we 
computed the cluster mass relation, rj, defined in that paper, 
but now included our exte nded list of cluste rs with X-ray 
counterparts. As defined in I Clark et al.l l|200ct ). r/ relates the 
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Figure 4. /-ba nd Antennae HST image showing positions of X- 
ray sources from lZezas et al.l ||2006| 1. North is up and east is to the 
left. Red squares are optical counterparts to X-ray sources with 
IR cluster counterparts identified in this paper. Black circles are 
those X-ray sources lacking cluster counterparts. 



function of X-ray detections per mass as a function of cluster 
mass and can be formalized with the following equation: 



Nx(M c ) = N C i(M c ) ■ V (M C ) ■ M c 



(1) 



The quantity, Nx{M c ) is the number of detected X- 
ray sources with an IR cluster counterpart, Nci(M c ) is the 
number of detected clusters, and r)(M c ) is the fraction of X- 
ray sources per unit mass, all as a function of cluster mass, 
M . 

If rj(M c ) increases or decreases over a range in M c , 
this means there could be something peculiar about mas- 
sive cluster physics to favour or suppress XRB formation. 
In contrast, a constant rj(M c ) across all M c would indicate 
that more massive clusters are more likely to have an XRB 
simply because they have more stars. 

We estimated cluster mass using 7\" s -band luminosity, 
Mk b , but expressed in terms of the observational quantity, 
flux, Fk s ■ Given a set of r\ values plotted against cluster 
flux, a significant slope could indicate an over-abundance of 
X-ray sources in more massive clusters, more than would 
be expected from simple scaling arguments. Fig. 6 shows 
T](Fk s ) computed for each cluster, but binned by bins of 
Fk s — 4.1 x 10 6 DN _1 . As can be seen, there does not 
appear to be a significant slope and rj is consistent with a 
constant value of ?](Fk s ) = 5.5 x 10 -8 with an uncertainty 
of q w = 7.0 x IP" 9 an d is in agreement with the value found 
in lClark et all (|2008T l. 

To estimate the strength of this result, we performed 
a x 2 test between the mean value o f rj and a fitted li ne to 
the plotted values of t) (Fig. 6) (see lClark et al.ll200St ). We 
found a AE;^ 2 = 0.45. Since any value of AEx 2 less than 
one indicates no significant slope, this suggests there is no 
slope in r\. This differs from the value of AE^ 2 = 1-9 when 
only using the 15 clusters found with associations to X-ray 
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Figure 5. (J — X s ) vs K s colour-magnitude diagram for all clus- 
ters in Antennae. Clusters with X-ray sources are designated by 
the luminosity class of the associated X-ray source. LLX: Lj < 
3 X 10 38 ergs s~\ HLX; 3 X 10 38 ergs s" 1 < L x < 1 X 10 39 ergs 
s _1 , and ULX: L x > 1 X 10 3 9 erg s -1 . Blue clusters with X-ray 
sources are from IClark et al.l l|2007l ) , while red clusters with X- 
ray sources are the additional sources added in this work. The 
dashed line signifies the photometric cut-off we set for statistical 
purposes. 



sources (|Clark et alj 2008) and seems to strengthen the case 
that there is only a fiat slope in the functional form of rj(M c ). 



3.2 Spectral Evolutionary Models 

During our spectrophotometric analysis of cluster proper- 
ties, we fit models to cluster colours across all six IR and 
optical bands. We estimated cluster colours by performing 
photometry using the same sized aperture across all bands. 
We defined a constant photometric aperture as ~3<r of the 
J-band Gaussian PSF, where the full width at half maxi- 
mum (FWHM) is 1.2 arcsec. The background annulus had 
a radius of ^6-10cr of the PSF. Considering HST resolved 
many of the IR sources into multiple components, this large 
photometric aperture encompassed these conglomerations. 
Since we expect that these sources are part of the same, 
larger structure, it is appropriate to include photometry of 
them in spectral evolutionary models. 

A bright , 2MASS star was used to compute J and K s 
magnitudes |Brandl et al.l 120051 ) . We derived HST magni- 
tudes using zero point s listed in table 28.1 of the HST Data 
Handbo ok (IVoitlll997l). Applyi ng colour transformations de- 
fined in iHolt zman et all (|l995f ), we converted all HST mag- 
nitudes to Johnson UBVI magnitudes. We expressed errors 
in magnitude, a m , as tjfi ux divided by the mean flux. In the 
case of the optical filters, a m consists of the additional er- 
rors in the zeropoint and colour transformations, all of which 
were added in quadrature. 

After performing photometry on all X-ray source coun- 
terparts seen across all IR and optical bands, we eliminated 



Figure 6. This figure displays r](Ffc g ) plotted versus M^. The 
bins are Fk s = 4.1 X 10 6 DN -1 in size. Error bars are the un- 
certainties in the mean value of r\ added in quadrature with the 
Poisson uncertainty in each bin. The dotted line is the mean of 
the five r\{Fj^ s ) values. 
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Figure 7. Example of a typical \ fit between a cluster counter- 
part and a BC model. We include the derived cluster properties 
and the quality of the fit parameters. Notice the small residuals 
between the model and the data in the lower graph. 
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Figure 8. Summary of results from \ 2 fitting to BC models. 
Error bars are ranges in cluster parameters for models within 
plus one of the best-fitting \ 2 sum. We divide these plots by X- 
ray source luminosity into three separate regions: LLX: Lx < 
3 X 10 38 ergs s~\ HLX: 3 X 10 38 ergs s" 1 < L x < 1 X 10 39 ergs 
s — 1 , and ULX: Lx > 1 x 10 39 ergs s _1 . Notice that there is no 
obvious trend between cluster properties and the luminosity of 
the associated X-ray source. 



six sources suffering from poor sky subtraction due to crowd- 
ing, thus resulting in negative flux measurements. There- 
fore, we continued further analysis on 21 counterparts to 
X-ray sources. Performin g a \ 2 minimization tech nique, we 
fit Bruzual-Charlot CBC; iBruzual fc Chariot] 120031 ') spectral 
evolutionary models to the cluster magnitudes available for 
each source to determine mass and age. We fixed the metal- 
licity to Z= 0.02, whi ch is an accepted value for young clus - 
ters in the Antennae (|Christoperll200Sl ; iBastian et "aTl l2009'). 
These parameters gave us a more complete understanding 
of the XRB environments in the Antennae. Selecting the 
model with the best-fitting \ 2 value, we selected all models 
that had a % 2 within one, % 2 + 1> of this value. The resulting 
fits gave us a range in age and reddening for each cluster. 
We estimated reddening by iteratively removing Ay from 
uncorrected magnitude cluster SEDs. By iteratively picking 
values for Av in the range 0.0 - 3.0 mag and in steps of 0.1 
mag, we selected the Ay that contributed to our best-fitting 
model. By allowing the value of Av to float, we accounted 
for measured uncertainties in reddening. Each fit is listed 
in Table 3 and an example fit is shown in Fig. 7. We only 
include those fits with Xn < 1-6, which is the la confidence 
interval in Xn f° r the number of degrees of freedom in our 
model fits, n — 5. Xn ls equal to x 2 / n - 

The shift in magnitude between the model and the data 
(A) gave us a mass estimate for each IR/optical counterpart. 
A contains information on the distance modulus and mass 
of the cluster. Subtracting off the distance modulus to the 
Antennae, rrid = 31.4 mag (for Ho = 75 km s _1 Mpc^ 1 ), 
we were left with a difference in luminosity between our 



cluster and a 1 Mq cluster as listed by the BC models. 
Converting this difference into a change in flux gave us the 
cluster mass in solar masses. Before computing this change 
in flux, we renormalized the luminosity difference, expressing 
it in terms of magnitudes instead of colours We renormalized 
the luminosity difference to Mk b for all fits, wh ere Mk, was 
comp uted using PSF, aperture photometry (see lClark et al.l 
l2007ri . 

We summarize the results from all model fits in Fig. 
8. We found ages between ~ 9 X 10 6 - 2 X 10 7 yr. Masses 
ranged from 7 x 10 5 - 2 x 10 e Mq, with most ~ 10 6 Mq. 
The extinction varied from Ay = 0.3 - 1.5 mag. Finally, we 
did not find a trend between these cluster properties and 
the luminosities of their associated X-ray sources. 



4 CONCLUSIONS 

In this work, we made a strometric frame-t ies between Chan- 
dra X-ray coordinates (|Zezas et al.ll2006h . WIRC IR pixel 
positions, and HST optical coordinates to search for star 
cluster counterparts to X-ray point sources. Us i ng th e list 
of 120 X-ray sources in table 3 of IZezas et all (2006), we 
found 38 of these sources are associated with an IR counter- 
part, 35 of which are confirmed star clusters. This expanded 
on previous work in which we found IR star cluster asso- 
ciations to 1 5 of th e 49 X-ray source s listed in table 1 of 
IZezas et al.l |2002T ) (|Clark et al.ll2007h . Our new sample of 
IR cluster counterparts includes 11 of the previousl y iden- 
tified 15 counterparts discussed in lClark et al.l l|2007h . Both 
this past study, and what we present here, indicate that most 
X-ray sources, roughly two thirds, do not spatially coincide 
with a cluster. A comparison between X-ray luminosity and 
the separation from the associated cluster, does not show a 
trend, thus indicating no relation between XRB type and 
where it forms in it's parent cluster. 

Extending our frame-tie from the X-ray to the optical 
using the IR as an intermediary, we were unable to find re- 
liable optical counterparts due to the crowded HST fields. 
This problem of identifying counterparts to X-ray sources 
was also encountered in a sim ilar study, which focused on 
optical counterparts to ULXs l|Ptak et al.| [2006'). Thus, in- 
stead, we found the optical counterparts to 27 of the IR star 
clusters associated with an X-ray source. An examination 
of the positions of the cluster counterparts in the Antennae 
indicates that the majority are in the spiral arms and bridge 
region between these two galaxies. This seems to clearly in- 
dicate that those X-ray sources with counterparts are tied 
to star formation in these interacti ng galaxies. 

Using photometry taken from lBrandl et all (|2005T ). we 
made an IR, JK 3 photometric study of the 33 IR star cluster 
counterparts whose X-ray sources have listed luminosities. 
We found most of the IR cluster counterparts are bright, ~16 
mag in K s and in slightly redder clusters, (J — K s ) = 1.1 
mag as compared to (J — K s ) — 0.8 mag for the general 
p opulation of cluste rs. This confirms our results presented 
in IClark et al l i|2007h , in which we performed a similar study 
with the smaller sample of IR counterparts. 

Following the work in I Clark et alj l|200ct ). we explored 
the relationship between c luster mass and the detected num- 
ber of X-ray sources. In IClark et all (|2008l ) we defined a 
function r\ relating the number of X-ray detections per mass 
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Table 3. Bruzual-Charlo t Model Fits. Chandra Src ID numbers follow the naming conven- 
tion of lZezas et alj (|2006|). Xn i s the normal \ 2 statistic divided by the degrees of freedom 
for the fit; five for all filters. Q is a goodness-of-fit measure, showing the probability that 
X 2 statistic exceeds a particular value of x by chance. 



Chandra Src ID Log(Ljf) 


Log(Age/yr) 


A 


V 


Log(M/M ) 


A. n 


Q 


11 


38.34 


6.900- 


-7.477 


0.94- 


-1.84 


5.37- 


-6.11 


0.71 


0.47 


27 


39.77 


7.000- 


-7.477 


0.98- 


-1.28 


5.05- 


-5.58 


0.62 


0.54 


38 


38.23 


7.040- 


-7.440 


1.04- 


-1.54 


5.91- 


-6.33 


0.94 


0.32 


51 


38.31 


6.900- 


-7.360 


0.48- 


-0.58 


5.98- 


-6.48 


0.44 


0.70 


83 


38.60 


6.740- 


-8.306 


0.70- 


-2.80 


5.88- 


-7.22 


0.18 


0.93 


87 


37.83 


6.940- 


-7.420 


0.71- 


-1.11 


5.63- 


-6.17 


0.60 


0.56 


89 


37.34 


7.000- 


-7.477 


1.48- 


-1.78 


5.13- 


-5.66 


0.58 


0.58 


95 


38.75 


7.220- 


-7.360 


0.71- 


-0.81 


6.17- 


-6.30 


1.47 


0.12 


101 


38.02 


7.120- 


-7.160 


0.45- 


-0.45 


5.97- 


-5.99 


1.10 


0.24 


102 


37.95 


6.920- 


-7.440 


0.35- 


-0.65 


4.51- 


-5.08 


1.39 


0.14 



as a function of cluster mass. We found r\ is consistent with 
a mean value of i](Fk s ) = 5.5 x 10 -8 . Using a \ 2 test, we 
compared a fitted slope of the plotted tj(Fk s ) values to the 
mean ti(Fk s ) value. We found a AEy 2 = 0.45. Consider- 
ing we found a AE X 2 = 1.9 in lClark et all 1)20081) . our new 
study shows a stronger relation between t)(Fk b ) and the 
mean value, effectively ruling out any inclination towards 
more X-ray sources residing in more massive clusters other 
than through simple scaling arguments. 

Including the 27 star cluster counterparts to X-ray 
sources seen across all UBVIJK a ba nds, we fit Bruzual- 
Charlot (BC; lBruzual fc Chariot] I2003T I spectrophotometric 
models to 10 of these clusters. The BC model fits indicate 
the X-ray-source-associated clusters are 7xl0 5 -2xl0 6 M 
in mass, ~ 9 x 10 6 - 2 x 10 7 yr in age, with extinction vary- 
ing between Av = 0.3 - 1.5 mag. These proprieties indicate 
star cluster counterparts to X-ray sources in the Antennae 
tend to be young and massive, which are consistent with 
th ese X-ray sources being associated with star formation. 
In lPtak et al.l (|2006l ). these authors also found similar mass 
results for their counterparts to ULXs. 

While we can use multiwavelength photometry to de- 
scribe cluster properties, there remains some uncertainty in 
these characteristics due to errors in magnitude and model 
limitations. In future work, we plan to acquire spectra of An- 
tennae cluster counterparts and refine estimations of their 
properties. 
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